The Random Field Ising Model on Hierarchical Lattices II: Ground State Critical 

Properties 
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The ground state critical properties of the Random Field Ising Model (RFIM) on the diamond hi- 
erarchical lattice are investigated via a combining method encompassing real space renormalization 
group and an exact recurrence procedure. The local magnetization and the nearest neighbors pair 
correlation function are exactly calculated. The fixed-point joint probability distribution of cou- 
plings and local fields are numerically obtained and analyzed, indicating that the critical behavior 
of the model is governed by the zero temperature disorder fixed point. The critical exponents asso- 
ciated with the order parameter and correlation length are estimated showing an universal behavior 
regarding the choice of the initial probability distribution for initial local fields being the symmetric 
continuous Gaussian or discrete delta-bimodal. 
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I. INTRODUCTION 

The universality of the critical behavior of the Ran- 
dom Field Ising Model (RFIM), with finite dimension, 
has been object of controversy. While some authors in- 
dicate that the critical behavior of the RFIM, in cubic 
lattices, does not dependppn the considered random field 
probabilities distributionEl, others argue that the RFIM 
with Gaussian or bimodal distributions belong to differ- 
ent universality classesa - such as it occurs for the. RFIM 
in infinite dimensions (mean field approximation p. Both 
the worksBH calculate exactly the ground state of the 
RFIM and apply the finite size scaling approach to get 
the critical exponents. According to the real space renc*- 
malization group analysis proposed by Bray and MooreQ, 
which argues that the critical behavior is governed by a 
disorder fixed point at T = 0, these results are valid for 
finite temperatures. However, very recently, Duxbury 
and MeinkeB have shown that, within the mean field ap- 
proach, the RFIM presents distinct behaviors between 
the ground state and the finite temperature case. If mean 
field approach conclusion is applicablc_to finite dimen- 
sion systems, the results of Swift et aln and Hartmann 
and Nowakcl will not reflect the RFIM behavior for finite 
temperature. However, the mean field approach can be 
pathological and not valid for finite dimension systems 
if, for instance, the local field fixed point distribution is 
GaussianQ. 

In this paper, which follows a previous study of the 
phase -diagram and thermodynamic potentials of the 
RFIMB (hereafter I), we study the ground state of the 
RFIM, either with the ±-ffo delta-bimodal probability 
distribution for the random fields or the Gaussian one, 
in diamond hierarchical lattices. We use a procedure 
that encompasses the real space renormalization group 
approach and an exact recursive process allowing us to 
get the local magnetizations and the nearest neighbors 
pair correlation functions. In the next section, we will 
describe the model and we present the renormalization 



equations properties. The fixed point distribution is ob- 
tained and analyzed in section III , where we comment on 



the applicability of the ground state results for systems 
with finite temperatures. In section IV, we calculate the 



ground state critical exponents for the RFIM using two 
distinct methods and compare their results. Finally, in 
section M, we summarize our results. 



II. THE MODEL 



The RFIM Hamiltonian is given by: 



H 






Z^HjCTi, 



(1) 



where J represents the ferromagnetic coupling between 
spins located at sites i and j of a diamond hierarchical lat- 
tice and Hi is the random magnetic field that acts on the 
Ising spin er,; = ±1. As usual, < ij > indicates that the 
sum must be done only on first neighbors. The quenched 
random fields variables can be chosen from a Gaussian 
or delta-bimodal probability distribution VH {Hi), with 
zero mean and variance Ho. 

The diamond hierarchical lattice with N generations, 
on which the model is defined, is constructed by replac- 
ing each bond of the N — 1 hierarchical lattice by the 
generating cell, as depicted in figure [l| The process be- 
gins with one bond joining the root sites. This lattice 
scaling factor equals two and its graph fractal dimension 
depends on the number of connections (p) of the gener- 
ating cell, that is dj = log(p)/log(2). For the p = 3 and 
p = 4 cases, considered in this paper, the dimension of 
the lattice is 2.58 . . . and 3, respectively. 

Applying the real space renormalization group ap- 
proach, the T — coupled renormalization equations for 
the coupling constants and fields are obtained by deci- 
mating the internal sites of the generating cell as sketched 
in figure ffl, giving: 




FIG. 1: The diamond hierarchical lattice is built in a recursive 
process (decoration). Starting from the primitive cell (left), 
we replace each bond of the preceding hierarchy by the basic 
cell (right). The decimation process is done in the spirit of 
the real space renormalization group in order to obtain the 
renormalization equations (J2a|) to (|2q). 




prise the ferromagneticJsing model with uniformt3t 
aperiodic interactions^], the spin-glass Ising mode, 
and the ferromagnetic Potts model defined on the di- 
amond hierarchical lattices, as well as the Ising model 
(both the pure and the spin-glass cases} defined on the 
Wheatstone bridge hierarchical latticaia and on the m- 
sheet Sierpinskii Gasket fractal latticed. 

For the RFIM ground state the coupled recursive equa- 
tions between local magnetization and pair correlation 
functions are written as 



< a > = A x + A 2 < Mi > 

+ A 3 < fi 2 > +A 4 < M1M2 > (3a) 

< cr^i > = A 2 + Ai < m > 

+ A 4 < fi 2 > +A 3 < fnn 2 > (3b) 

< <7[I 2 > = A 3 + A A < m > 

+ Ax < fi 2 > +A 2 < muz > (3c) 



where, 



1 b 

-max(|ifj|, |Jij - J 2i \)] 

1 b 
H' = H + -J2[^M\Ju\,\Hi + J 2i \) 

-max(|Jii|, \Hi - J 2l \)] 

1 b 
H' = H + -^[max(|J 2i |,|i? i + J u \) 



i=l 



-max(|J 2l |, \Hi - Jx, 



(2a) 



(2b) 



(2c) 



Contrary to the finite temperature case (to cite our 
work), in T — the renormalization equations are linear. 
In this way, we can divide all the bonds and local fields 
by a common factor at each renormalization step, with- 
out modifying the physical properties of the model. This 
property is extremely useful when analyzing the renor- 
malization flow for two main reasons: firstly, it avoids 
some numerical overflows and secondly choosing the nor- 
malization factor as the width of the renormalized field 
distribution, we keep this distribution with fixed width 
(unitary) and we get the fixed point distribution. Thk 
approach was adopted previously by Cao and MachtaEl 
for another three dimensional hierarchical lattice. 

The method for evaluation of the local magnetization, 
which is summarized in I and fully presented inl3, is a gen- 
eralization of an approach developed to study the local 
magnetization of the Ising spin glass model on diamond 
hierarchical lattices under an external magnetic fieldD. 
This approach has also been extensively applied to inves- 
tigate the local order parameter of several magnetic mod- 
els defined on hierarchical lattices. Those studies com- 



AAx = sgn(h + Kx + K 2 ) + sgn(h + K X - K 2 ) 

+ sgn(/i - Kx+ K 2 ) + sgn(h - K x - K 2 ), 

AA 2 = sgn(h + Kx + K 2 ) + sgn(/i + K x - K 2 ) 

- sgn(/i -K x + K 2 ) - sgn(h -K^- K 2 ), 
AA 3 = sgn{h + Kx + K 2 )- sgn(h + K x - K 2 ) 

+ sgn(h - Kx+ K 2 ) - sgn(/i - K x - K 2 ), 
4A 4 = sgn(h + Kx + K 2 )- sgn(/i + K X -K 2 ) 

- sgn(h - Kx+ K 2 ) + sgn(h - K x - K 2 ). 



(4a) 

(4b) 

(4c) 

(4d) 
-1 if 



In equations (^), sgn(x) = 1 if x > 0, sgn(x) 
x < and sgn(x) = if x = 0. 

The equations ([|) can be used to calculate the whole 
set of local magnetizations and nearest neighbors pair 
correlation functions, as a function of the strength of the 
random fields, by means of an inflation process, as il- 
lustrated in figure |]. From these data we can calculate 
directly certain thermodynamic potentials, such as the 
average magnetization and the internal energy. 



III. FIXED POINT DISTRIBUTION 

Under the renormalization process, the initial distri- 
bution of local fields and the single delta function distri- 
bution for ferromagnetic couplings evolves either to the 
ferromagnetic fixed point, if Hq < Hq, or to the param- 
agnetic one if Hq > Hq. In an appropriated distribution 
parameter space, the flow of the renormalized distribu- 
tions towards these "fixed-points" defines the respective 
basins of attraction. Within the frontiers of these basins 
of attraction, one should find the unstable fixed point 
distribution, which governs the system critical behavior. 
Therefore, the determination of the fixed point distribu- 
tion is an important issue to study the critical behavior 



of the model. In particular, the form of the fixed point 
distribution of the random fields determines whether the 
criticality is actually dictated by the ground state fixed 
point (as argued by the renormalization group )u or the 
result of the mean field approach remains valid in finite 
dimensions. 

To analyze the renormalization flow, we generate a 
pool of 10 5 triplets {Jij,Hi,Hj} from the initial dis- 
tributions, each one corresponding to the coupling con- 
stant and the local fields acting of a given bond (see I). 
We use equations (0) to generate a renormalized pool of 
triplets with same size. The renormalization sequence 
of triplets pools determines the renormalization flow. In 
the paramagnetic phase (Ho > Hq), the couplings av- 
erage between spins vanishes while the width of the lo- 
cal field distribution (a H =< H 2 > 1 / 2 ) grows, but with 

< J > /ch —* 0. On the other hand, in the ferromag- 
netic phase (H < Hq), the coupling average grows faster 
than the local fields, so that < J > /an — » oo. In the 
critical point (Ho = Hq), however, the transition between 
the two assintotic behaviors occurs, that is < J > /an 
it is unstablely constant. Therefore, to determine the 
fixed point distribution we must find the initial width of 
the random field probability distribution for which the 
renormalized < J > /an remains constant. Using this 
procedure, we found the critical fields for the three di- 
mensional lattice (Hq/J = 0.904 ±0, 001 for the bimodal 
distribution and H§/J = 0.997 ± 0,001 for the Gaus- 
sian one) while, for the lattice with dimension 2, 58 ... , 
we obtain H§/J = 0.587 ± 0, 001 for the bimodal distri- 
bution and H§/J = 0.610 ± 0, 001 for the Gaussian one. 
These latter values confirm that, at least for the diamond 
hierarchical lattices, the ferromagnetic long range order 
is stable for dimensions bellow d = 3. Furthermore, we 
have also confirmed that, for d = 2, there is no phase 
transition (even an infinitesimal random field distribu- 
tion width drives the renormalization flow toward the 
paramagnetic phase). 

To generate the fixed point distribution, one may keep 
the distribution in the fixed point applying, at each step 
of the renormalization, a small perturbation to the bond 
distribution, diminishing about 20% the difference of 

< J > of the renormalized distribution to its value at 
the fixed point distributions. In figure g we show the fixed 
point distribution when the initial random field distri- 
bution is the bimodal one. A very similar plot is ob- 
tained for the Gaussian initial distribution. It is worth 
to notice that the fixed point distribution is correlated. 
This correlation corroborates the necessity to consider 
the triplets {Jij,Hi,Hj} instead of uncorrelated fields 
and bonds distributions. In contrast of the result got- 
ten by Cao and Machtala, the fixed point distribution for 
the diamond hierarchical lattice, the bond distribution 
does not present finite probability in the zero limit of the 
bond strength, however (see figure ||) it also presents an 
asymmetry around < J >. Within this aspect, coij re- 
sults are similar to that obtained by Newman et a/Ji3. In 
the figures and 0, we present the integrated fixed point 



distribution for the coupling constants and the random 
fields, respectively, the latter beina,very close to a GauSj 
sian distribution. Cao and MachtaEl and Newman et alx3 
also obtained fields distributions very close to the Gaus- 
sian one. Therefore, despite the bond distribution de- 
pendence on the hierarchical lattice, the fixed point field 
distribution seems to be rather universal. 
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FIG. 2: Integrated fixed point probability distribution for 
the delta-bimodal initial random field distribution with width 
Hg/J = 0, 903. 



To test numerically the Gaussian shape of the (inte- 
grated) fielcL-distribution, we perform the Lilliefors test 
of normalityta. Intuitively, this test measures the " dis- 
tance " between the accumulated probability distribu- 
tion obtained from the renormalization process and that 
of the Gaussian. The result of the test shows that the 
fixed point distribution for the random fields is, at least, 
very close to a Gaussian, considering a significance level 
of 0.1E3. Therefore, accordingly to referenced, the results 
of the mean field approximation are not applicable to 
the finite dimension system and can we can determine 
the critical behavior of the RFIM by studying its ground 
state. 



IV. CRITICAL EXPONENTS 

To calculate the critical exponents of the RFIM we 
use two distinct methods. The first one is based on the 
way that the joint probability distribution V(H, J) moves 
away from the fixed point distribution through the renor- 
malization process, while second one is a finite size scaling 
analysis of the magnetization. 
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FIG. 3: Integrated fixed point bond probability distribution 
for the delta-bimodal initial random field distribution with 
width HS/J = 0, 903. 



FIG. 4: Fixed point local field probability distribution for 
the delta-bimodal initial random field distribution with width 
HS/J = 0, 903. 



Following Bray and Mooraa, we define the exponents x 
(associated with the renormalization of an infinitesimal 
symmetry breaking field), y (associated with the growth 
of the width of the distribution) and z = 1/v (y being 
critical exponent associated with the correlation length). 
These three exponents can be directly calculated from tbp 
fixed point distribution, as described in the referencea. 
Thus, we can calculate 



*(&i) 



log 6 



log A 
log 6 



and 



_ _ 1 _ log(h'/h) 
v log 6 

where the primes are relate the renormalized quantities, 
[. . . ] is a configurational average, b = 2 is the lattice 
scaling factor, A is the average of the rescaling factor of 
the distribution (that is, an) and h = uh / J — {?H / 'J) > 
(ajj/J)* being the value of <th/J at the fixed point. Us- 
ing these equations, we get the exponents presented in 
table 1. The exponents a and (3 were calculated through 
the hiperscale relations proposed by Bray and MooreQ. 

As shown in table |, all the critical exponents calcu- 
lated here, either for three dimensional case or to the 
d = 2, 58 . . . one, indicates that the critical exponents do 
not depend on the choice of the initial probability dis- 
tribution. We emphasize that the critical exponent p is 



Exponent 


Bim 

d = 2.58. . . 


odal 

d = 3 


Gaussian 
d = 2.58. . . d= 3 


X 


2.59 ±0.02 


2.99 ±0.02 


2.62 ± 0.04 


2.98 ± 0.02 


V 


t.292 ± 0.003 


1.489 ± 0.002 


1.293 ± 0.003 


1.488 ± 0.002 


V 


3.41 ±0.06 


1.996 ± 0.008 


3.41 ±0.04 


2.00 ±0.01 


a 


-2.41 ±0.09 


-1.02 ±0.02 


-2.41 ± 0.06 


-1.02 ± 0.02 


p 


0.0 ±0.2 


0.02 ±0.04 


-0.1 ± 0.2 


0.04 ± 0.04 



TABLE I: The critical exponents x, y e v were calculated by 
the analysis or the renormalization flow. p| e /? were calculated 
with the hiperscale relations proposed in 



very small, being possibly zero, for both probability dis- 
tribution cases. This almost zero value of /3 may indicate 
a discontinuity in the magnetization. The exponent a, 
by it turn, is negative, showing that the specific heat do 
not to diverge. These results agree with the-maiority of 
the previous renormalization group resultsaOEj. How- 
ever, our results present an appreciable difference in the 
value of the exponent v when compared with, the value 
v = 1 obtained by Boechat and ContinentinocJ. The ori- 
gin of this discrepancy is related with the uncorrelated 
distribution of bonds and local fields considered by these 
authors. Therefore, we verify once again the importance 
of the correlation in the critical properties of the system. 
To calculate the exponents v and f3 through finite size 
scaling, we use the equations (||) to calculate local mag- 
netizations and, then, we calculate the average magneti- 



zation 



1 



N a 



1 Ns 

I = 1 



< <r<"> > 



(5) 



where N a is the number of samples and N s the number 
of sites of the lattices. The absolute value of the magne- 
tization of a single sample is taken because the random 
field breaks the symmetry in each sample arbitrarily to 
the positive or negative magnetization state. Therefore, 
the average of the module of the magnetization is the 
relevant parameter. 

We get the average magnetization for different lattice 
sizes (3 to 8 hierarchies), as shown in the figured. Tak- 
ing a linear regression (represented by the straight lines 
in the figure 0), for each lattice size, in the region where 
the magnetization begins to fall abruptly until the begin- 
ning of the finite size effects, we find the point where this 
curve intercept the x-axis. We define this point as the 
"critical field " for that lattice size. Although arbitrary, 
this criterion is rather robust, since the inclusion or ex- 
clusion of some points in limits of the regression region 
does not modify very much the value of the exponents. 
On the other hand, the finite size scaling predicts that 
all the pseudo-critical fields must converge to the critical 
field with the same critical exponent iMB, 




where N it is the lattice size. 

Therefore, plotting log [(H$(N) - H%/ H§] x log(iV) 
(figure @), we estimate the critical exponent v as the 
slope of the best fit of the plot. With the value of v 
on hands, we obtain the value of /3 from the scaling law 
M ~ N-P/ v f((H Q - H^M 1 /^). In the table || we 
present our results for v and (3 obtained by both con- 
sidered methods and observe that, within the error bars, 
the critical exponents are equal. 
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FIG. 6: Illustration of the calculation of the u critical expo- 
nent. The circles refers to the Gaussian distribution, while 
the diamonds label the bimodal distribution results. 



Approach 


Exponent 


Bimodal 


Gaussian 


Renormalization 
Flow 


V 


1.996 ±0.008 


2.00 ±0.01 





0.02 ± 0.04 


0.04 ± 0.04 


Finite size 
scaling 


V 


2.0 ±0.1 


1.9 ±0.1 





0.02 


0.02 



TABLE II: Comparison between the critical exponents ob- 
tained via the renormalization flow and finite size scaling. 



FIG. 5: Ground state magnetization for the RFIM with Gaus- 
sian random field probability distribution. The lines are the 
linear regression used to calculate the critical field for a given 
lattice size. From above to bellow the curves correspond to 
lattice with 3 to 8 hierarchies. 



(Hg{N)-HS)/H5~N- 1 ' u , 



(6) 



V. CONCLUSION 

In this paper, we study the RFIM on diamond hi- 
erarchical lattices, using the real space renormalization 
group approach. Studying the renormalization flow, we 
obtained the joint fixed point probability distribution 
for coupling and fields. The analysis of the local field 
fixed point distribution shows that this distribution is 
(at least) very close to the Gaussian one. Therefore, the 



critical properties of the model are dictated by the dis- 
order critical point (ground state), An contrast to what 
happens in the mean field approach^. 

We calculate the critical exponents of the ground state, 
showing that the RFIM with delta-bimodal and Gaussian 
distribution belong to the same universality class, either 
for the three dimensional lattice or the d = 2,58... lat- 
tice. We found an almost zero (3 critical exponent, in- 
dicating a discontinuity in the magnetization. A scaling 
analysis with larger lattices could elucidate this question. 
However, such analysis exceeds our computational capac- 
ities in the present moment. 



Acknowledgments 



We thank to M. Contincntino and J.F. Fontanari 
for helpful discussions. AR is greatful to CNPq and 
FACEPE (Brazilian granting agencies) for financial sup- 
port. This work also received financial support from 
FINEP (under the grant PRONEX 76.97.1004.00), CNPq 
and CAPES. 



1 M. R. Swift, A. J. Bray, A. Maritan, M. Cieplak, J. R. 
Banavar, Europhys. Lett. 38 273 (1997). 

2 A. K. Hartmann and U. Nowak, Eur. Phys. J. B 7 105 
(1999). 

3 A. Aharony, Phys. Rev. B 18, 3318 (1978). 

4 A. J. Bray and M. A. Moore, J. P hys. C 18 L927 (19 85). 

5 P. M. Duxbury and J. H. Meinke, |cond-mat/0012042J (To 
appear in Phys. Rev. E). 

6 A. Rosas and S. Coutinho, cond-mat/0107404 (submited 
to Phys. Rev. B). 

7 O. D. da Silva-Neto, DSc. Thesis, Universidade Federal de 
Pernambuco, (1999). 

8 M. S. Cao and J. Machta, Phys. Rev. B 48, 3177 (1993). 

9 M. E. Fisher, Phys. Rev. 113, 969 (1951). 

10 W. A. M. Morgado, S. Coutinho and E. M. F. Curado, J. 
Stat. Phys. 61, 913 (1990). 

11 W. A. M. Morgado, E. M. F. Curado and S. Coutinho, 
Rev. Bras. Fis. 21, 247 (1991). 

12 E. Nogueira J r., R. F. S. Andrade and S. Coutinho, cond 



n at 0008194 (to appear in European J. Phys B (2001)). 



E. Nogueira Jr., S. Coutinho, F. D. Nobre, E. M. F. Curado 

and J. R. L. de Almeida, Phys. Rev. E 55, 3934 (1997). 

E. Nogueira Jr., S. Coutinho, F. D. Nobre and E. M. F. 

Curado, Physica A 257, 365 (1998). 

G. Camelo-Neto, MSc. Thesis, Universidade Federal de 

Pernambuco, (1999). 

R. de P. A. Lima, MSc. Thesis, Universidade Federal de 

Pernambuco, (1999). 

M. E. J. Newman, B. W. Roberts, G. T. Barkema and J. 

P. Sethna, Phys. Rev. B 48 16533 (1993). 

E. J. Dudewicz and S. N. Mishra, Modern Mathematical 

Statistics (John Wiley & Sons, 1988). 

A. Falicov and A. N. Berker and S. R. McKay, Phys. Rev. 
B 51 8266 (1995). 

B. Boechat and M. A. Continentino J. Phys. 2 5277 (1990). 
M. N. Barber, Finite size scaling in: Phase Transitions 
and Critical Phenomena, 8, edited by C. Domb and J. L. 
Lebowitz (Academic Press, New York, 1983). 



